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Abstract. We present our model to interpret the SED of galaxies. The model for 
the UV to sub- mm SED is already well established (Silva et al., 1998). We remind 
here its main features and show some applications. Recently we have extended the 
model to the radio range (Bressan et al., 2001), and we have started to include the 
X-ray emission from the stellar component. 



1. Modelling the UV to sub- mm SED of galaxies 

We have developed a model (GRASIL, http://grana.pd.astro.it) for the 
SED of galaxies including dust effects (Silva et al., 1998; Silva, 1999). 
We consider a geometry more complex and realistic than previous 
treatments, and take into account 3 distinct environments where stellar 
radiation interacts with dust (Figure 1): (1) molecular clouds (MCs) 
+ young stars (clumped component), (2) diffuse ISM (cirrus) + 
older stars (diffuse component), (3) AGB envelopes (included in the 
SSP library, Bressan et al., 1998). Therefore both stars and dust have a 
realistic 2-phase distribution. The association between young stars and 
MCs is modelled to mimic the birth of stars in MCs and their gradual 
escape as they grow old: the fraction of each SSP energy radiated inside 
MCs is decreased with its age (Figure 1). It is the first time that such 
an age-dependent or selective extinction is modelled (see Granato et 
al., 2000 for interesting applications). 

The dust model consists of normal grains in thermal equilibrium, 
small grains fluctuating in temperature and PAH molecules. 



2. Modelling radio emission in star-forming galcLxies 

We have extended our model to the radio. The thermal radio emission 
is proportional to the production rate of ionizing photons: 
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The non thermal radio emission is the sum of a dominating term, 
proportional to the supernova rate (cahbrated on the Galactic value) 
and one due to supernova remnants (based on the E— D relation): 



erg/s/Hz Vl.4y ' Vl.4 

See Bressan (these proceedings) and Bressan et al. (2001) for details 



3. Applications 

3.1. Fitting the SEDs 

As an example, the fit to the starburst galaxy M82 is shown in Figure 2. 
Note that the inclusion of the radio emission potentially (with the 
adopted calibration) helps in constraining the starburst age, due to the 
different timescales involved in the IR (ongoing SFR, last few Myrs) 
and radio (recent SFR, a few 10^ yrs ago) emission. 

3.2. The redshift of sub-mm selected galaxies 

The sub-mm to radio spectral index has been recently used to estimate 
the redshift of sub-mm selected galaxies (e.g. Carilli & Yun 2000). 
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Figure 8. Left: Fit to the SED of M82. The starburst age and e-folding time are 
respectively 25 Myr and 10 Myr (continuous line) and both 50 Myr (dotted line). 
An older starburst is compatible with the UV to sub-mm SED but not with the 
radio emission. A free-free absorption 4QHz ~ '^'^ been applied. Right: Fit to 
the SED of Arp220. 




Figure 3. Left: 850/im-1.4 GHz spectral index vs z for the ARP220 model observed 
at different ages during the starburst. Right: linear combination of the (850/im-1.4 
GHz) and (8.4-1.4 GHz) spectral indexes vs z for the same models. 



Exploiting our model it is possible to investigate the dependence of this 
index on the evolutionary status of galaxies. In Figure 3 the 850/im- 
1.4 GHz spectral index is plotted vs redshift using the model fitting 
the SED of ARP220 (Figure 2, this seems to be representative of the 
SCUBA galaxies, e.g. Barger et al., 1999), but for different ages of 
the starburt. The redshift dispersion around z=2 (the typical value 
estimated for SCUBA galaxies) is Az = 1.8. To reduce this dispersion, 
we propose to use a radio-radio spectral index (e.g. 8.4-1.4 GHz), almost 
independent of z, in combination with the previous one. Indeed, with 
the linear combination shown Figure 3, the Az due to the evolution 
becomes 0.5 around z=2. 
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Figure 4- Left: Predicted evolution with age of the contribution of different stellar 
components to the X-ray emission of SSPs. The minimum mass for BH formation 
is 14M0 and Po = 0.02 s. Right: Predicted stellar X-ray SED for M82 based on the 
UV-to-radio fitting model. Data from Moran & Lehnert (1997), Cappi et al. (1999). 

4. Stellar X-ray emission from galaxies 

Recently we have started to include in our model the stellar contribu- 
tion to the X-ray emission of galaxies. Our starting point has been the 
work by Van Bever & Vanbeveren (2000). They find that the dominat- 
ing stellar components in the X-ray for starburst galaxies are massive 
binaries with a black hole or a neutron star as a primary, and an OB 
star as a secondary either in the H-shell burning phase or on the main 
sequence, and pulsars. We have included these components, and the 
contribution from supernova remnants, in the SSP library. The main 
uncertainties in the predicted output are (Figure 4): the binary fraction, 
the minimum mass for BH formation and the initial spin period Pq of 
pulsars (whose Lx oc Pq'^)- A further uncertainty is the SED to assigne 
to each component. A first attempt for the stellar X-ray SED expected 
from the UV-to-radio fitting model for M82, is shown in Figure 4. 
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